NASA Technical Memorandum 102320 


A Thermoelastic Transversely Isotropic 
Thick Walled Cylinder/Disk 
Application; An Analytical 
Solution and Study 


S.M. Arnold 
Lewis Research Center 
Cleveland, Ohio 


October 1989 


(NASA-TM-102320) A THERHOEL AST IC 
TRANSVERSELY ISOTROPIC THICK WALLED 
CYLINDER/DISK APPLICATION: AN ANALYTICAL 
SOLUTION AND STUDY (NASA) 46 p CSCL 20K 

G 3/39 


N90-12950 


Unci as 
0240384 





A THERMOELASTIC TRANSVERSELY ISOTROPIC THICK WALLED CYLINDER/ 
DISK APPLICATION: AN ANALYTICAL SOLUTION AND STUDY 


S.M. Arnold 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135 


ABSTRACT 


A continuum theory is utilized to represent the thermoelastic behavior of a thick 
walled composite cylinder that can be idealized as transversely isotropic. A multi axial 
statement of the constitutive theory employed is presented, as well as, the out of the plane 
of isotropy, plane stress and plane strain reductions. The derived analytical solution 
presented is valid for a cylindrical tube or thin disk with a concentric hole, subjected to 
internal and/or external pressure and a general radial temperature distribution. A specific 
problem examined is that of a thick walled cylinder subjected to an internal and external 
pressure loading and a linear radial temperature distribution. The results are expressed in 
nondimensional form and the effects on the response behavior are examined for various 
material properties, fiber orientation and types of loadings. 


INTRODUCTION 

Recently, a continuum theory was presented for representing the thermoelastic 
behavior of unidirectional composites which can be idealized as a locally transverse 
isotropic homogeneous continua [1]. Applicability of a continuum theory for predicting 
deformations of composite (textured) materials depends relatively upon characteristic 
structural dimensions, the severity of gradients (i.e. stress, temperature, etc.), and the 
relative scale and periodicity of the material’s internal structure. 

Here, application of this transversely isotropic theory is made to the plane strain 
problem of a thick walled cylindrical tube subjected to internal and external pressure and a 
radial temperature distribution. The corresponding problem of a thin disk (plane stress 
case) with a concentric hole is also considered. Motivation for selecting the thick walled 
cylinder problem for investigation was twofold. The first is that it provides a simple 
structure that can be solved analytically, thereby supplying a benchmark problem for 
structural analysis code verification. Secondly, the thrust chambers of reusable rocket 
engines, notably the NASA Space Shuttle main engine, and tubular composite test 
specimens may be geometrically idealized as such. This enables a qualitative examination 
of the effect of fiber reinforcement on the stress, strain and deformation distributions. 
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The present study will be confined to examining the effect of material properties, 
fiber orientation and type of loading (in a general non-dimensional setting) on the 
thermoelastic behavior of a cylinder so as to develop an intuition with regard to the effects 
of transverse isotropy. In particular four types of pressure loadings (i.e. interior only, 
exterior less than interior, uniform and exterior greater than interior) and three types of 
linear radial temperature profiles (i.e. inner surface hotter than outer, uniform , and outer 
surface hotter than inner) will be examined. 

This study begins with a brief description of the thermoelastic transversely isotropic 
constitutive theory and its plane strain and plane stress simplification. Next the analytical 
solution is presented in general form and then in a normalized form for a specific linear 
radial temperature distribution. A limited parametric study is then presented. 


CONSTITUTIVE EQUATIONS 

Given a linear hyperelastic material reinforced by a single family of fibers the stress 
and strain components can be related through a normality structure utilizing a 
complementary energy function Q( cr. j,d. dj,T), i.e. 


dn 
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Here cr- j denotes the components of (Cauchy) stress, c-j the components of infinitesimal 
strain and the symmetric tensor d-dj is formed by a self product of the unit vector dj, 
which denotes the local fiber direction and T represents temperature. 

As discussed in [1] the complementary energy function may be expressed as: 
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Note that in this report, the complementary energy function of equation (2) has been 
extended to include the effect of temperature. 

Substituting equations (2) — (4) into equations (1) and associating coefficients 
A,B,C,D,E,F and H with physically measurable parameters (e.g. Young’s modulus and 
Poisson ratio) as specified in [1] results in the following linear thermoelastic stress strain 
relations. 


£ ij 




E™ E 


ij 


* {3( £ ~I L)I 1 + + irs_ ]I 4 * <¥j 

-T) * 4i 


T "L 
* ( 1 

— = 1 V 


IJ 


2G t 


ij 


(5) 


where 


and 


'ij 




— denotes the Kronecker delta function 
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G^, — Shear modulus (transverse) in the plane of isotropy 

— Shear modulus (longitudinal) in a plane normal to the plane of isotropy. 

Erp — Young’s modulus (transverse) in the plane of isotropy. 

E^ — Young’s modulus (longitudinal) normal to the plane of isotropy. 

^rp — Poisson’s ratio (transverse) that characterizes the transverse strain reduction 
in the plane of isotropy due to a tensile stress in the same plane. 

— Poisson’s ratio (longitudinal) that characterizes the transverse strain 

reduction in the plane of isotropy due to a tensile stress in a direction normal 
to it. 

— characterizes the thermal expansion normal to the plane of isotropy 
Orp — characterizes the thermal expansion in the plane of isotropy 
AT = T 2 — Tj ; TpT 2 are two distinct temperatures. 

TWO DIMENSIONAL SIMPLIFICATION 

Simplification of the multiaxial linear thermoelastic stress — strain relations of (5), 
is made to that of plane stress and plane strain. Out of the plane of isotropy (see figure 1), 
the plane stress assumption requires that the stress traction 

r i = Vj = 0 



while plane strain requires 


£ U = 0 (7) 

provided the local fiber direction (d-) and nj (indicating the directionality of the plane 
stress or strain assumption) are defined as, 
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( 8 ) 


dj = cos ip = 1 

d 2 = sin tp n 2 = 0 

d 3 = 0 n 3 = 0 

Note that the shear strain e ^ is i n general considered to be nonzero, while would have 
been required to be equal to zero, if d^ had been nonzero. Given the above assumptions 
the following expressions relating the nonzero stress and strain components are obtained. 


Plane Stress: 

fj. = (7 r /Erji + {— J'rp/Erp +[l4p/E,p — I^/Ejjs 2 }^ 

+ OrpAT (9) 
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Plane Strain: 


e r = cr r {l/E T — C 2 / A} + {- v^f E x +[i/ t /E t - vj Ejs 2 

— CB/A}«T0 + {oj, -CD/ A} AT (13) 

e & = {1/E t + [1/2 G l — l/E T -i/ L /E L ] S 2 + 

[1/E t +(1+2*/ l )/E l -1/G l ]s4 — B 2 /A (14) 

-f" { I^rji/Efji -f (l^rji/Erj-i — i^/Ej^)S 2 — BC/A )(^4' 

[arp+ (a^-Qrp)s 2 — BD/A ]AT 

e z6 = {[1/2G L -1/E rt /E L ](1-B/A)cs + 

[1/B t +(1+2i/ l )/E l -1/G l ](c8>-(B/A)cSs)}^ 

+ {{u t /E t -u l /E l )cs -C/A([1/2G l -1/E t -i/ l /E l ]cs 
+ [l/E T +(l+2i/ L )/E L -l/G L ] C 3s)}«7 r 
+ [(c^-a^cs -(D/A){[l/2G L -l/E T -i/ L /E L ]cs 
+[1/E t +(1+2i/ l )/E l -1/G l ]c3s)}]AT 

(15) 

* 2 =-(B/A)^-(C/A)cr r -(D/A)AT (16) 

where A,B,C and D are defined in Appendix A and 
s 11 — sin n y> 

n = 1,2,3 and 4 

n n 
c = cos tp. 

ANALYTICAL SOLUTION 

The problem under consideration is the solution of a thick walled circular cylinder 
with end caps composed of a locally transversely isotropic material (unidirectional 
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composite) which is subjected to a uniform internal (Pj) and/or external (P Q ) pressure and 
a radial temperature distribution, see fig 2. The inner and outer radii are denoted by a and 
b , respectively, while T & and denote the applied temperature along the inner and outer 
radii of the cylinder. Generally a cylinder is classified as being thick walled when the wall 
thickness exceeds the inner radius by more than approximately 10%, i.e. b— a > 0.1 or A = 
b\a > 1.1. 


General Solution 

The governing equations for this axisymmetrical problem, are: 
Equilibrium — neglecting body forces 



Strain — Displacements: 


d u 
dT 


6/i = 


u 


and 


Compatibility: 
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where 


<7., <7-0 is the radial and tangential (hoop) stress 
respectively 

e r ,60 is the radial and tangential (hoop) strain 
respectively 

Substituting the appropriate stress-strain relations , equations (9-10) and (13-14) 
into equation (19), simplifying and applying equilibrium the following Euler— Cauchy 
equation, describing the radial stress, results: 

=/d(r*K)) -Q(r)) 

dr dr 
or 

r 2 d 2 a r + 3r d a i + (1 - 1 ) a = Q(r) (20) 

TFT TF J x 

Here /3 and Q(r) are dependent upon the constitutive equations employed, e.g. plane stress 
or plane strain assumption (see Appendix A). Implied in equation (20) is the assumption 
that the material parameters (i.e. Ep,Ep,yp,i/p,Gp,a:p and Op) are not strongly 
temperature dependent and therefore can be taken to be independent of radial location. 
Inclusion of this temperature dependence would only contribute higher order effects in the 
solution and distract attention from the primary emphasis of examining the effect of 
anisotropic coefficients of thermal expansion. 

The general solution of (20) is 

a r =C 1 r m i + C 2 r m 2 + Q 1 (r)-Q 2 (r) (21) 
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where 



rm ' / r“( m .+ X )Q(r) dr 
m 1 — m 2 


Q 2 (r) = r m 2 / r“< m ! + J >Q(r) dr 

m 1 — m 2 


m 1 = — 1 + E 


m 2 = — 1 — E 


E = 



and as is always positive, the roots (m^ and m 2 ) are therefore real and distinct. The 
general solution for the tangential stress is obtained by substituting equation (21) into 
( 17 ); 


< j q = 1) r 1 + C 2 (m 2 + 1) r 2 + (m^ 1) Q x (r) 

- (m 2 + 1) Q 2 (r) (22) 


The constants and C 2 are found by applying the boundary conditions: 


K>r=a = - P i 
K>r=b = ~ P o 
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to equation (21). The resulting expressions are 
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where S(r) = Q^r) - Q 2 (r). 


Substituting and C 2 into equations (21) and (22) yield the following expressions 
representing a ■ and Cq for either the plane stress or plane strain assumption and an 
arbitrary radial temperature distribution. 


= [ a - ” 2 Pi - IT” 2 P 0 ] r m 1+ [b“” 2 a" ll ' 1 " 2 P„ - a’” 2 b">“ m2 P ) J r m 2i 


(b' 


mi- m 2 _ m i*~ m 2 

a 


(b 1 


m I— m 2 _m i“ in 2 

cl 


) 


{a Hn2 Q 1 (a)-b Hn2 Q 1 (b)}r ini + {a mi_m2 b Hn2 Q 1 (b)-a“ m2 b nil_ra2 q 1 (a)}r 


m2 


(b 


m m 2 mi- m 2 

a- 


) 


(b' 


mi - m 2 


- a 


m [— m2 


) 


+ Qi(r) 


{a“ ra 2 q 2 (a)-b^ , 2 q 2 (b)}r rai _{a raiHn2 b Hn 2 q 2 (b)-a Hn 2 b nliHn 2 q 2 (a)}r in2 


(b’ 


m 1- m2 


- a 


mi- m 2 


) 


(b' 


mi — m2 


- a 


m i“ m2 


) 


w 


10 


( 24 ) 
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If a linear radial temperature distribution is assumed 


AT = T -T + (T b“ T a>(r-a) 

1 0 

and the resulting radial and tangential stress are non-dimensionalized, then the following 
expressions are obtained: 
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Note II, Z Q and Z ^ are given in Appendix B for the plane stress and plane strain case. 
Neglecting all temperature dependence, the above equations can be shown to agree with the 
work of Lekhnitskii [2] 

Knowing the non-dimensional radial and tangential stress components, the 
non-dimensional radial deformation, shear strain and either the longitudinal strain or 
stress components can be determined utilizing equations (6), (7), (18), (26) and (27). 

The resulting expressions are: 
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where 



and 


T = T, /T 
b' a 

T = T /T 
o' a 

and ^,S,C , ) £>,Z2,Zg ) Z^,Zg ) Zg,Z^,Zg,Zg,and Z^q are defined in Appendix B. 


Isotropic Simplificatior 


Under isotropic conditions (i.e. E=E^=Erp, G=G^=G,p v=v^=Ur^ and a= a^=arp) 
the above plane strain expressions simplify to the following, which agree with solutions 
found in Timoshenko [2]. 
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'zO= 0 


(35) 


(ff r ) 4 ("p 


_ (Ea AT) 


(36) 


and P,A,p,T and AT are defined as before. 


DISCUSSION AND RESULTS 

Comparing, in analytical form, the transversely isotropic solution, equations (26) 
—(31), to their isotropic counterparts equations (32) —(36), the following general 
observations can be made as to the significant role this special anisotropy plays. 

1) The shear strain (eq. 29) is present for an arbitrary angle <p, however 
when the material is isotropic or the preferred direction is along the axis of 

symmetry (<p=0°) or perpendicular (ip=90 o ) no shear strain is computed. 

2 

2) The radial and tangential stress distribution is at most a function of (1/p ) 

2 

and A for an isotropic material, whereas with an anisotropic material 
dependence on p and A is now a function of material properties and 
angle of orientation. 

3) In an isotropic material the longitudinal stress (cr ) is independent of radial 
location (p) and equal to a constant when only mechanical loads are present 

Conversely when the material is anisotropic a is dependent upon material 
properties, angle of orientation and p. A similar statement can be made 
about the longitudinal strain (e ) when the plane stress assumption is 
invoked. 
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4) Considering only thermal loading, no stress is developed under a uniform 

temperature field if the material is isotropic. However stresses are developed 
in an anisotropic material provided the angle of orientation is not parallel 

with the axis of symmetry (^=0°). 

Note similar statements are applicable concerning items 1,2 and 4 when the plane stress 
assumption is invoked. 

To gain a better understanding of the role anisotropy plays with regard to the stress 
distribution; various loading conditions, material parameters and fiber orientations will be 
considered. The presentation will be divided into three parts, i) internal/external pressure 
loading, ii) thermal loading and iii) thermomechanical loading. All results are presented in 
non-dimensional form. 


Pressure Loading: 

Considering only the application of internal and external pressure, equations (26) 
and (27) reduce to the following: 


[1-PA 1+S ] ^E-l + [PA 1+S - A 2S ] ^-(1+E) 
(A 25) -l) A 22 - 1 

53 [1~P^ 1+S ] _ S [P^ 1+S - ^ 2S ] -(UZ) 

0 2S -1) Jl 2S - 1 


(37) 


(38) 


Taking the derivative of the normalized radial and tangential stress with respect to 
the normalized radius ( p ) and equating them to zero allows one to obtain critical points 
were the stress distribution would remain stationary. The critical location for the radial 
stress is 
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( 39 ) 


1 

'(1+S) (PA 1+X -A 2S )]^S 

(E-1)(1-PA 1+S ) 


and for the tangential stress 


P 


cr 


1 

■_(1 + E)(PA 1+S -^ 2S ')]^E 
(S-l) (1-PA 1+S ) 


(40) 


For p cj . to be within the cylinder or disk it must be positive valued and be within 
the limits 1< p^< A. Clearly no realistic stationary point is available for the tangential 
stress. A realistic stationary point may, however exist for the radial stress distribution as 
indicated in (39) depending upon the material parameters, fiber orientation, thickness and 
pressure ratio. 

Calculations were performed to illustrate the above dependency on a transversely 
isotropic material with the assumed base material parameters listed in Table 1, for four 
types of pressure loadings, i.e.; 

1) Interior pressure only, P= 0 

2) Exterior pressure less than interior, Pel 

3) Equal interior and exterior pressure, P=1 

4) Exterior pressure greater than interior, P>1 

in figures 3 through 6, respectively. These figures show the effect fiber orientation, i.e. <p= 
0,30,60 and 90°, has on the normalized radial stress distribution for a thick walled cylinder 
(plane strain assumption), with ratio of outer to inner wall radius (A) of 5.0. 

Observe that when tp=0°, the fiber orientation is aligned parallel with the axis of 
symmetry or equivalently the material is isotropic, i.e. E=l, thus we find from (39) that 
P cr = od, therefore no stationary points exist within the cylinder or disk. Thus the 
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maximum/minimum radial stress occurs either at the inner or outer surface, see equation 
(37); 


if p = 1 

and 

if p = A 



depending upon the values of the applied pressures. This is consistent with well known 
classical solutions [3,4]. 

For example in figures 3 (P=0) and 4 (P<1) the maximum stress occurs at the inner 
surface, while in figure 5 (P=l) it is constant and in figure 6 (P>1) the maximum stress 
occurs along the outer surface. Clearly the development of stationary points within the 
cylinder is dependent upon orientation and loading. Case in point is observed in figure 4, 
were for <p=6 0 and 90°the radial stress actually becomes less than the applied exterior 
pressure, which utilizing isotropic intuition should be the minimum stress obtainable. 
Conversely, the maximum stress continues to remain either at the interior or exterior 
surface independent of orientation, as one might expect since the stress component is 
governed by boundary conditions, i.e. the type of loading. The corresponding tangential 
stress distribution is shown in figures 7-10 respectively, for P=0, P<1, P=1 and P>1 and 
as indicated by (40) no stationary points are observed regardless of fiber orientation. 

Work with isotropic thick walled cylinders [4] has shown that when both internal 
and external pressure act simultaneously the location of the maximum tangential stress, be 
it along the inner or outer surfaces, is dependent upon both geometry and the ratio of 
external to internal pressure. This conclusion is supported by the present work, figures 
7—10, with the additional dependency of material property and fiber orientation due to the 
anisotropy of the material. This is evident in the following equation which expresses the 
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ratio of tangential stress at the inner (p=l) and outer (p=A) surfaces: 


S = 


a 


a 


_inner = 

0 , 
outer 


( 1+A 2S ) - 2PA 1+S ) 
2 A S_1 - P ( 1+A 2S ) 


( 41 ) 


In Table 2, equation (41) is utilized to obtain the required pressure ratios which would 
produce the following four example tangential stress ratios; 

1) uniform tangential stress (S=l) 

2) equal and opposite (S=-l) 

3) inner tangential stress zero, outer tangential stress maximum (S=0) 

4) 0<S<1, which insures maximum tangential stress at the outer surface and 
that the sign of Ug remains the same throughout. 

The fact that the sign of the tangential stress, as well as other stress components, can 
change from tension to compression or vice versa is important and should be considered 
when the elastic or inelastic material behavior is different in tension and compression, e.g. 
concrete or ceramics, as premature failure may result. 

An example of the variation of the longitudinal stress (a ) with radius depending 
upon fiber orientation is illustrated in figure 11 for the case of P=0.5. The associated 
normalized radial displacement is shown in figure 12. Figure 13 shows the distribution of 
shear strain, e^g, with angle orientation for the case when P=2.0. Note that for <p = 0 and 
90° the shear strain is zero. 


Variation of Material Property: 

Thus far the effect of loading and fiber orientation on the distribution of stress and 
displacement has been examined. The above equations also indicate, however that the 
strength of anisotropy (specific material properties) will affect this distribution and as the 
material properties in Table 1 were chosen arbitrarily, the effect of their variation will now 
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be examined. As indicated in [1] to ensure convexity of the complementary energy function 
certain restrictions must be imposed on the material parameters. These restrictions are 
restated in non-dimensional form in Table 3. Next the effect of changing k and when 
t'rj, is held fixed at 0.3 is investigated. 

Figures 14 through 19 illustrate the effect that changing material properties have on 
the radial and tangential stress distributions for the case when the exterior pressure is 
twice that of the interior (i.e. P=2.0) and the fibers are oriented at <p=6 0°. In each figure 
only one material parameter is varied, while the remaining are equal to those prescribed in 
Table 1. Figures 14 — 19 clearly indicate that varying k has the most significant effect 
upon both the radial and tangential distribution ( with regard to curvature and 
maximum/minimum values) while varying £ has the least. This trend appears to be 
independent of the type of pressure loading. 

Thermal Loading: 


Now considering only thermal loading equations (26) and (27) reduce to 


! _ {A 1+E -l} E-l + {A 1+S -A 253 } — (1+E) 


{A 2S - 1} 


{A 2S - 1} 


, _ + {^ 2£ } ,-( 1+E) 


{A 2S - 1} 


{A 2E - 1} 


(42) 


1 - p^ 1 - > ,-( 1+£ ) 


{A 2E - 1} 


{A 2S - 1} 


2, _ Et A 2+E -l} ,8-1 _i{A 2+S -A 2S } ,-(l + I) 
” 2E 


{A 2S - 1} 


{A - 1} 
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( 43 ) 


Again taking the first derivative of the normalized radial and tangential stress with 
respect to the normalized radius (p) and equating them to zero provides a condition to 
determine the critical location for the radial stress, i.e. 


7 (E - 1) \n jl+S 7 ,2 + E frj 7 \"| 53 — 2 

l-V + V - ( Z o + Z l)J P 

(r L -i) 

- [z 0 a 1+s + z^ 2+s - (z 0 + z^a 25 ] P _E_2 = o 


(A 251 -!) 


and tangential stress 


(44) 


2Z 


S(E—l) rjj il+S n ]2 + E /rj rj \ 1 

i ~ J LV + V ~ ( z o + z i)J p 


53-2 


(A 2E -1) 


[Z 0 ^ 1+S + Z 1 A 2+S — (Z q + ZjJA 251 ] ^ 2-2 = 0 

(A 22 -l) 


(45) 


where Z Q and Zj are defined in Appendix B, for the plane stress or plane strain 
assumptions. A general solution for p critical is not straight forward, although it is clear 
that p cr depends upon geometry (A), material properties (e.g. £, k, and a) , angle 
orientation (ip) and thermal load (T and T). Under isotropic conditions (i.e. 53=1, a=l) 
and assuming plane strain, equations (44) and (45) yield 


2V|_ (46) 

A+l 
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P = 


-A z 

~X+1 


(47) 


respectively. Clearly no stationary point exists for the tangential stress, as p would be 
imaginary. The radial stress may, however have a stationary point depending of course 
upon the ratio A. 

Calculations were performed to investigate the above dependency for the 
transversely isotropic base material of Table 1, when 

1) inner surface is hotter than the outer surface (T<1) 

2) uniform temperature field (T=l) 

3) outer surface is hotter than the inner surface (T>1). 

For simplicity the temperature at the inner surface (T ) is assumed to be equal to the 
reference temperature (T q ), i.e. T = 1 , except for case 2 (T=l) as both Zq and would 
be equal to zero and no temperature loading would exist if T = 1. 

Figures 20 —23 show the radial and tangential stress distributions through the 
thickness of the cylinder at 30° increments in fiber orientations, for thermal loadings of T 
< 1 and T > 1 respectively. Clearly the sign for both the maximum radial and tangential 
stress (be it tensile or compressive) is dictated by the type of loading (i.e. T <1 or T>1) 
for a fixed T. Also the location of the maximum radial stress is now angle and material 
dependent since the radial stress at the inner and outer surface are zero. The maximum 
tangential stress occurs at the inner surface independent of angle and thermal loading. 

The variation with radius for a z and e * are given in figures 24 and 25 for load case 
1 (T<1). As one would expect a varies linearly with p at a 0° fiber orientation and the 
normalized shear strain (e^) is zero. A non— zero shear strain and a nonlinear distribution 
of a is observed for off-axis orientations and (as with pressure loadings) the magnitude is 
highly dependent upon material properties, thickness (A) and angle. 
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The radial stress distribution produced for a given fiber orientation under a uniform 
temperature (T=l) distribution is shown in figure 26. No radial (or tangential) stress is 
calculated for a 0° orientation as one would expect since both < 7 . and Oq are in the plane of 
isotropy. Note that the sign (tensile or compressive) of the stress (radial or tangential) is 
dependent upon whether the applied thermal load is above (T<1) or below (T>1) the 
reference temperature. In figure 26, T is equal to 0.5. It is important to realize that both 
expansion and contraction could take place within the thickness of the cylinder if both T 
and T are less than or greater than one, but not equal. 

Variation of Material Properties: 

The effect of varying material property parameters, when the cylinder is subjected 
to a thermal load is very similar to that shown for the pressure loadings. For example in 
figures 27 and 28 the effect on the radial stress distribution of changing £, (i.e. £=2,5 and 
10) and k (i.e. k = 0.8, 1.0, 2.0 and 4.0) are shown respectively for the case when T=0.5, T 
= 1.0, a=0.5 and at a fiber angle of 60°, see figure 20. In addition to varying these 
mechanical material parameters, the coefficient of thermal expansion will play a significant 
role in two ways. 

The first is its effect on the sign of the stress; for example if a=a^/ a^> 1 (cf. figure 
26 where a=2.0 and T<1) is tensile whereas if a<l (cf. figure 29 where a=0.5 and T<1 
and T=l) a is compressive. The second is in the actual magnitude and how it affects the 
magnitude of the various stress components as illustrated in figures 30 and 31 where (J r and 

A A A 

a q versus p are displayed respectively. For the case when T=2.0, T=1 , <p= 60° and a = 
—0.25,0.0,0.25,0.5,1.0 and 2.0. Figure 30 and 31 clearly indicate that as a approaches zero 
a minimum radial and tangential stress state is reached. This a value, which produces a 
minimum stress state is not unique as it is highly fiber orientation dependent. Examples of 
this are shown in figures 32 and 33 for the identical case shown in figure 30 but with the 
angle <p being changed to 30 and 90° respectively. 
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It is important to realize that in a macroscopic sense the magnitude of a=a^/ap 
represents the ratio between the longitudinal and transverse coefficient of expansion or in 
an averaged sense the amount of a mismatch there is between fiber and matrix assuming 
longitudinally aligned fibers. For example assume we have a cube of composite material 
subjected to a change in temperature. If the amount of elongation parallel and transverse 

^ A 

to the fibers is the same then 0 ^= 0 ^ and a — 1 . Similarly an a = 1 would be equivalent 
to no a mismatch between the fiber and matrix provided that the coefficients of thermal 
expansion in all directions are equal [5]. The point is that elimination of any a mismatch 
between fiber and matrix may not be the most desirable path to follow if one wishes to 
minimize the macro stress state within a structure. Therefore by selecting the appropriate 
system one might be able to optimize the a mismatch for a particular structure, applied 
load(s) and fiber orientation. Clearly the implied assumption here is that the micro 
stresses developed due to this optimum a mismatch are insufficient to cause self 
destruction of the material. 

Thermomechanical: 

Thermomechanical results can be obtained from the above results by superimposing 
the various pressure loadings on those of the pure thermal loading. Equations 26 — 31 are 
valid for any arbitrary pressure and/or linear temperature loading. The response behavior 
is highly dependent upon the given type of loading, material properties, fiber orientation 
and geometry. An example is shown in figures 34 and 35 were the normalized radial and 
tangential stress versus the normalized radius p are shown respectively when a uniform 
pressure loading (cf figures 5 and 9 ) is superimposed on a temperature profile whereby the 
outer surface is twice that of the inner (cf figures 22 and 23). One observes from figure 34 
that the radial stress now changes sign from compressive to tensile stress and that the 
stress distribution is now relatively insensitive to fiber orientation. The insensitivity to 
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orientation is merely a coincidence as clearly shown in figure 36 where the temperature at 
the outer surface is now four times that of the inner. 

CONCLUSIONS: 

A continuum theory representing the thermoelastic behavior of unidirectional 
composites has been applied to the problem of a thick walled cylindrical tube subjected to 
internal and external pressure and a radial temperature distribution. The resulting 
analytical solution was then specialized for a linear radial temperature distribution and was 
expressed in non-dimensional form. A limited parametric study was conducted in which 
four types of pressure loadings and three types of radial temperature profiles were 
examined to develop intuition with regard to the effects of transverse isotropy. 

This analytical solution provides a benchmark problem for structural analysis code 
verification and may be employed to illustrate the stress, strain and deformation 
distributions within a unidirectional tubular composite test specimen or the thrust 
chambers of reusable rocket engines. The results clearly illustrate that the stress and strain 
response is highly dependent upon geometry, fiber orientation, material properties and type 
of loading. Finally it was illustrated that to minimize the macro stress or strain state 
within this structure subjected to the given pressure and thermal loadings a non unique a 
mismatch would be desired. This mismatch is highly dependent upon fiber orientation, 
material property and geometry. 
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APPENDIX A 

Material dependent parameters of the Euler— Cauchy Equation. 


Plane Stress (utilizing equations (6)) 


P = W - + ( - -- ~^ l )s 2 + (1 + ( 1+2t/ i,) _1 ) s 4 } 


e t g l e t e l 


e t e l g l 


Q(r) = E T{ (a L -o T )s 2 AT + + (fl L -a T )s 2 ) r 1_ AT } 

/? dr 


Plane strain ( utilizing equations (7)) 


/»- { * ♦(!-?- ♦ ( l ♦ (^i) . [ >s 4 - 1 } 

e t g l e t e l e t e l g l { 

[i - t 2 ] 

e t i 

Q(r) = (j’r^ AT 

dr 

where 


, -[(» L -» T )s 2 + (C - B)D ] 

P 1 = A 

{ A * ( 1 - ^ - £^L)s 2 + (1* W -h s 4 -B 2 } 
e t g l e t e l e t e l g l ^ 
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s BD j 

r 


{ 1 + ( 1 _2 _ 2 i / L)s 2 + ( 1 + ( 1 + 2 i / l ) _ 


®L E T e l 


Em E t 




and 


a = { i +(l-!-!!i)c 2 + (l + ( ^i ) -l)c 4 } 


*1 G l e t e l 


*T E L C L 


= { - "L t ( 1 + ( 1+2 "L) - 1 )s 2 c 2 } 

h h % h 


= {-!t + (!t--!i)c 2 } 


E| Erp Ej^ 


D 

AT 


— £ fl/j a T )c J 

= I(r) -I 0 


• n 


s = sm <p 


c n = cos 11 p 


n = 2,4 
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APPENDIX B 

Non-dimensional material dependent coefficients 


With the following defined quantities: 


E t 
« = T 


a= a L 


i = T o 


Em 


T = T b 


{ 1 + («-2-^L)c 2 + (1 — k + ) c 4 } 

{ -h + ( 1 + -«)s 2 c 2 } 


{ -"T + ( u l-^l)c 2 } 


= [ 1 + (a -l)c 2 ] 


sin 11 ip 


c n = cos 11 <p 


n = 1,2, 3, 4 


the non-dimensional material dependent coefficients for plane stress and strain assumption 
may be obtained. 
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Plane stress (utilizing equations (6)) 


S 2 = 1/0 

0 = { 1 + ( ft — 2 — 2 ^L)s 2 + ( 1 + ( 1+2t/ I,) — *)s 4 ) 

Zq = _E T a T T a [ (a - l)s 2 ] [ (A-T) — T ] 

PiO? - i ) 

Z 1 = ~ E T a T T a [1 + 2(a — l)s 2 ] (T-l) 

PiW " 1 ) wv 

Z 2 = (^ t ~^ l )s 2 -i/ t 

r 

Z 3 = 0 

Z 4 =l + (o-l)s 2 
Z 5 =(^ T -^ L )cs 

r 

Zg = ( E — 1 — ^L) cs + ( 1 + ( ^ +2y I.) — «) cs 3 

* t t 

Z^ = (a-1) cs 

Zg = + (Vf -^L)c 2 

Z 9 = ~ v l + ( 1 + f 1 + 2t/ L) - *) s 2 c 2 
Z 10 = 1 + ( a-1 ) c 2 


30 


Plane strain (utilizing equations (7)) 


E 2 = [1-^] jfi 

P = { 1 + ( k-2- 2v I)s 2 + ( 1 + ( 1+2i/ l) -k)s 4 -^} 

7 _T i 


Z Q = E T a T T a ^ [ (a - l)s 2 + {C-B)D] [ (A-T) -T] 

p i (^ - i + f) ~~i i^ry 

i 

Z 1 = -E T a T T a [1 + 2(n - l)s 2 + (C-2B)D] (T-l) 

p ± (4^ - i + f) ~r xx=ry 

i 

7 / \ 2 BC 

Z 2 " (*'T”*'l) s ~ l 

T 

h - P 


Z 4 = 1 + (a— l)s 2 — BD 

T 


h = (^ T -^r)cs-f{(*-l-^L)cs+ (1 + ( 1+2t/ L) -*)c 3 s } 


Z 6 = ( K - 1 - f/ L)cs (1 - B ) + ( 1 + ( 1+2v l) - *) (cs 3 - Vs) 
6 2 J 1 —f 1 

Z 7 = (a-1) cs-^{( E -l- *L)cs + ( 1 + ( 1+2 V -«)c 3 s } 
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TABLE 1 

Base Material Parameters 


E t = 15000 ksi 

" L = 08 
= 0.3 

a T = 3-OxlO -6 


£ = 5.0 
k = 2.0 

a= 0.5 


TABLE 2 

Required Pressure Ratio to Produce the Circumferential Stress Ratio 


S = -l 


S = 0 


S = 1 


P = 


P = 


i 1 2E 01 
1+1 +21 

~ ,2E 

1 + 1 +21 

l + l 22 
21 1+2 


E-l 

T+E 


P = 


21 2-1 - 1 - l 22 
1 + l 22 - 21 1+2 


1 + l 22 < p < 21 2-1 - 1 - l 22 

21 1+S 


0 < S < 1 


1 + l 22 - 21 1+2 


TABLE 3 

Inequalities Required to Maintain Convexity of Potential 

f >0 : ;>**: t > 2(/ L : *> 2 ( U,, L> 

L 1^, ' f 


( 

(i J'jK > 

, 2 

*> 2 ( 1 -l) 

2 

L 

■ £ ■ 

2.0 

0.836>i/ l 


*>1.8 

5.0 

1 .322>^ 


A e>0.72 

10.0 

1.87 >v^ 


/c>0.36 
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NOMENCLATURE 


— Complementary Energy Function 

— Vector denoting local fiber direction 

— total infinitestimal strain tensor 

— Cauchy stress tensor 

— Temperature 

— Physical Invariants k=l,2,3,4,5 

— Invariants k=l,2,4,5 

— Transverse shear modulus 

— Longitudinal shear modulus 

— Transverse Young’s modulus 

— Longitudinal Young’s modulus 

— Transverse poisson’s ratio 

— Longitudinal Poisson’s ratio 

— Longitudinal thermal expansion 

— Transverse thermal expansion 

— normalized thermal expansion coefficient, a^ctr^ 

— Stress traction vector 

— Angle between axes 1 and 2 or z and 0. 

— Internal pressure 

— External pressure 

— Normalized pressure, i.e P^/P Q . 

— inner wall radius 

— outer wall radius 

— temperature at inner wall 

— temperature at outer wall 

— ratio of temperature of outer wall to that of inner. 

— ratio of outer wall radius to that of inner. 

— radial displacement 

— normalized radial displacement, =u/a 

— radial location 

— normalized radial location, = r/a 

radial stress and strain component respectively, 
tangential stress and strain component respectively. 



Vz" 


r 

a e 


[z6 

T 


longitudinal stress and strain component respectively 

— normalized radial stress component. 

— normalized tangential stress component. 

— normalized longitudinal stress component. 

— normalized longitudinal strain component. 

— normalized shear strain component. 

— normalized reference temperature, = T /T 

Q oL 



Figure 2. - Schematic of thick walled cylindrical composite 
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Floure 3. - Normalized radial stress versus normalized radius 
for fiber orientations of 0, 30, 60 and 90° and interior pressure 
only, P* 0. 



Figure 4. - Normalized radial stress versus normalized radius for 
fiber orientations of 0, 30, 60 and 90° and interior pressure 
greater than exterior, P - 0.5. 




Figure 5. - Normalized radial stress versus normalized radius 
for fiber orientations of 0, 30, 60 and 90* and uniform 
pressure, P- 1. 


Figure 6. - Normalized radial stress versus normalized radius 
for fiber orientations of 0, 30, 60 and 90® and exterior 
pressure greater than Interior, P - 2 . 
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Figure 7. - Normalized tangential stress versus normalized 
radius for fiber orientations of 0, 30, 60 and 90° and 
interior pressure only, p - 0. 
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Figure 8. - Normalized tangential stress versus normalized 
radius for fiber orientations of 0, 30, 60 and 90° and 
interior pressure greater than exterior, P - 0.5. 



P 

Figure 9. - Normalized tangential stress versus normalized 
radius for fiber orientations of 0, 30, 60 and 90° and 
uniform pressure, P - 1. 


P 

Figure 10. - Normalized tangential stress versus normalized 
radius for fiber orientations of 0, 30, 60 and 90° and 
exterior pressure greater than Interior, P - 2. 
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Figure 11.- Normalized longitudinal stress versus normalized 
radius for fiber orientations of 0, 30, 60 and 90° and 
Interior pressure greater than exterior, P - 0,5. 



Figure 12. - Normalized radial displacement versus normalized 
radius for fiber orientations of 0, 30, 60 and 90° and interior 
pressure greater than exterior, P « 0.5. 



Figure 13. - Normalized shear strain versus normalized 
radius for fiber orientations of 0, 30, 60 and 90® and 
exterior pressure greater than interior, P » 2. 
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Figure 14. - Normalized radial stress versus normalized radius 
for varying EJ, i.e., 2, 5, and 10 at fiber orientation of 60® and 
exterior pressure greater than Interior, P « 2. 



Figure 15. - Normalized tangential stress versus normalized 
radius for varying £, I.e., 2, 5, and 10 at fiber orientation of 
60° and exterior pressure greater than Interior, P - 2. 



Figure 16. * Normalized radial stress versus normalized radius 
for varying k, I.e., 0.8, 1.0, 2.0 and 3.0, at fiber orientation of 
60® and exterior pressure greater than interior, P * 2. 
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Figure 17. * Normalized tangential stress versus normalized 
radius for varying ic, l.e.,0.8, 1 .0, 2.0 and 3.0, at fiber 
orientation of 60° and exterior pressure greater than 
Interior, P - 2. 


P 

Figure 18. - Normalized radial stress versus normalized radius 
for varying l.e., 1 .2, 0.8, 0.4 and 0.25, at fiber orientation 
of 60° and exterior pressure greater than interior, P » 2. 



Figure 19. * Normalized tangential stress versus normalized Figure 20. - Normalized radial stress versus normalized radius 

radius for varying u j_, l.e., 1 .2, 0.8, 0.4 and 0.25, at fiber ter fiber orientations of 0, 30, 60 and 90 s when the inner 

orientation of 60° and exterior pressure greater than interior, surface is hotter than the outer , T - 0.5. 

P * 2. 
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Figure 25. - Normalized shear strain versus normalized radius 
for fiber orientations of 0, 30, 60 and 90° when inner surface 
is hotter than the outer, T * 0.5. 
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Figure 27. - Normalized radial stress versus normalized radius 
for varying i.e. 2, 5, and 1 0 at fiber orientation of 60° when 
the inner surface Is hotter than the outer, T = 0.5. 
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Figure 26. - Normalized radial stress versus normalized radius 
for fiber orientations of 0, 30, 60 and 90° for a uniform 
temperature distribution, T ■ 1 .0, T « 0.5. 



Figure 28. * Normalized radial stress versus normalized radius 
for varying k, I.e., 0.8, 1 .0, 2.0 and 4.0, at fiber orientation of 
60° when the inner surface is hotter than than the outer, 

T * 0.5. 
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Figure 29. - Normalized radial stress versus normalized radius 
for fiber orientations of 0, 30, 60 and 90° for a uniform 
temperature distribution, T - 1 .0, with T * 0.5 and a - 0.5 
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Figure 31 . - Normalized tangential stress versus normalized 
radius for varying a, l.e., -0.25, 0, 0.25, 0.50, 1,0 and 2.0 at 
a fiber orientation of 60° when the outer surface Is hotter 
than the inner, T - 2.0. 
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Figure 30. - Normalized radial stress versus normalized radius 
for varying ft, l.e. ,-0.25, 0, 0.25, 0.50, 1 .0 and 2.0, at a fiber 
orientation of 60° when the outer surface Is hotter than the 
inner, T - 2.0. 
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Figure 32. - Normalized radial stress versus normalized radius 
for varying a, l.e., -0.25, 0, 0.25, 0.50, 1 .0 and 2.0, ata fiber 
orientation of 30° when the outer surface Is hotter than the 
Inner, T - 2.0 
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Figure 33. ■ Normalized radial stress versus normalized radius 
for varying &, l.e., -0.25, 0, 0.25, 0.50, 1.0, and 2.0 at a fiber 
orientation of 90° when the outer surface is hotter than the 
inner. T s 2.0. 


Figure 34. - Normalized radial stress versus normalized radius 
for fiber orientations of 0, 30, 60, and 90° when the outer 
surface is hotter than the inner, T * 2.0 and a uniform pressure, 
P-1.0. 


Figure 35. - Normalized tangentJaJ stress versus normalized radius 
for fiber orientations of 0, 30, 60, and 90° when the outer 
surface is hotter than the Inner, T - 2.0 and a uniform pressure, 
P-1.0. 


Figure 36. - Normalized radial stress versus normalized radius 
for fiber orientations of 0, 30, 60, and 90° when the outer 
surface is hotter than the inner, T - 4.0 and a uniform pressure, 
P- 1.0. 
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